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Abstract — The  robust  symmetrical  number  system  (RSNS)  is 
a  number  theoretic  transform  based  on  N  >  2  sequences 
that  can  extract  the  maximum  amount  of  information  from 
symmetrical  folding  waveforms.  The  sequences,  based  on  coprime 
moduli,  exhibit  an  integer  Gray  code  property  making  the  RSNS 
well-suited  for  many  applications  that  benefit  from  an  inherent 
error  detection  and  correction  capability  such  as  analog-to-digital 
converters,  direction  finding  arrays  and  radar  waveform  design. 
To  use  the  RSNS,  it  is  necessary  to  know  the  greatest  length 
of  combined  sequences  without  ambiguities,  called  the  dynamic 
range  M ,  for  which  only  a  few  closed-form  expressions  currently 
exist.  In  this  paper,  an  efficient  algorithm  for  computing  M  and 
its  position  within  the  combined  set  of  sequences  is  presented 
and  shown  to  be  independent  of  the  size  of  the  moduli.  The 
algorithm  is  used  to  generate  the  equations  for  several  groups 
of  additional  moduli  arrangements.  Closed  form  expressions  for 
M  are  conjectured  and  proved  using  the  obtained  congruence 
equations  that  define  the  ambiguity  locations. 


I.  Introduction 

THE  most  common  type  of  waveform  in  engineering 
science  is  the  symmetrical  folding  waveform  (e.g.,  si¬ 
nusoids).  Symmetrical  folding  waveforms  appear  naturally  in 
many  engineering  disciplines  and  system  analysis  techniques. 
To  extract  the  maximum  amount  of  information  from  symmet¬ 
rical  folding  waveforms,  symmetrical  number  systems,  each 
consisting  of  IV  >  2  integer  sequences,  were  formulated 
based  on  coprime  modular  systems.  Symmetrical  number 
systems  include  the  symmetrical  number  system  (SNS),  the 
optimum  symmetrical  number  system  (OSNS)  and  the  robust 
symmetrical  number  system  (RSNS).  To  effectively  utilize 
symmetrical  number  systems,  the  dynamic  range,  M,  defined 
as  the  greatest  length  of  distinct,  paired  sequences  (TV— tuples) 
that  contain  no  ambiguities,  as  well  as  its  beginning  position 
within  the  combined  A— sequences  must  be  known.  Closed- 
form  expressions  for  M  have  been  reported  for  the  SNS  and 
the  OSNS  in  [1]  and  [2],  respectively.  Unlike  the  SNS  and 
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OSNS,  a  general  closed-form  expression  for  the  RSNS  M 
does  not  exist  and  closed-form  expressions  have  only  been 
reported  for  a  limited  number  of  specific  cases  [3] — [5].  As  a 
result,  the  use  of  an  algorithm  is  required  to  determine  M  and 
the  beginning  position  of  the  sequence  for  moduli  sets  that  are 
not  covered  by  the  limited  number  of  cases. 

The  RSNS  has  an  inherent  integer  Gray  code  property 
that  makes  the  RSNS  particularly  attractive  for  error  control 
in  both  analog  and  digital  signal  processing  applications. 
The  RSNS  has  been  shown  to  be  useful  in  software  radio 
systems  for  sample  rate  conversion  [6],  and  in  electronic  [7], 
[8],  photonic  [9]  and  superconducting  [10]  folding  analog-to- 
digital  converters.  Due  to  the  inherent  symmetry  within  the 
modulus,  a  new  theoretic  transform  for  error  detection  and 
control  was  reported  in  [11]  and  applied  to  code  division 
multiple  access  wireless  communications  [12].  The  complexity 
of  direction  finding  antenna  systems  is  also  reduced  through 
use  of  the  RSNS  by  decomposing  the  spatial  filtering  operation 
into  a  number  of  parallel  sub-operations  [13].  Consequently, 
each  sub-operation  only  requires  a  complexity  in  accordance 
with  that  modulus  and  a  much  higher  spatial  resolution  is 
achieved  after  the  results  of  these  less  complex  sub-operations 
are  recombined.  The  use  of  the  RSNS  in  radar  waveform 
design  has  also  been  reported  in  [14]  to  extend  the  capabilities 
for  target  detection. 

In  this  paper,  we  present  an  efficient  algorithm  to  compute 
the  M  and  its  beginning  position  within  the  sequence,  for  a 
general  set  of  N  coprime  moduli.  The  algorithm  is  derived  by 
considering  the  location  and  distance  between  all  of  the  vector 
ambiguity  pairs  for  the  combined  A— sequences.  To  simplify 
our  derivation,  we  define  the  center  of  ambiguity  (CO A)  as 
the  midpoint  between  the  ambiguity  pairs.  Analysis  of  all  the 
ambiguity  pairs  leads  to  an  upper  bound  on  M  that  is  used 
to  improve  the  algorithm’s  efficiency.  An  A  =  3  example  is 
provided  to  demonstrate  the  steps  of  the  algorithm.  Also,  the 
algorithm’s  complexity  is  computed  and  compared  to  that  of 
a  naive  search  approach.  We  demonstrate  that  our  algorithm’s 
complexity  is  independent  of  the  moduli  size  and  represents  an 
improvement  by  several  orders  of  magnitude  when  compared 
to  the  naive  search  approach.  Our  algorithm  is  then  applied  to 
several  different  groups  of  moduli  sets  and  is  used  to  generate 
additional  closed-form  expressions  for  M.  The  closed-form 
expressions  for  M  are  developed  for  several  additional  A  =  3 
and  A  =  4  cases  and  are  verified  by  deriving  the  closed-form 


U.S.  Government  work  not  protected  by  U.S.  copyright. 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

2014 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-2014  to  00-00-2014 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


4.  TITLE  AND  SUBTITLE 

Extended  Closed-form  Expressions  for  the  Robust  Symmetrical  Number 
System  Dynamic  Range  and  an  Efficient  Algorithm  for  its  Computation 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

Naval  Postgraduate  School, Department  of  Electrical  and  Computer  report  number 

Engineering, Monterey, CA, 93943 

9.  SPONSORING/MONITORING  AGENCY  NAME(S )  AND  ADDRESS(ES )  10.  SPONSOR/MONITOR' S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

IEEE  Transactions  on  Information  Theory  60:3  (2014),  1-11. 

14.  ABSTRACT 

The  robust  symmetrical  number  system  (RSNS)  is  a  number  theoretic  transform  based  on  N  2  sequences 
that  can  extract  the  maximum  amount  of  information  from  symmetrical  folding  waveforms.  The 
sequences,  based  on  coprime  moduli,  exhibit  an  integer  Gray  code  property  making  the  RSNS  well-suited 
for  many  applications  that  benefit  from  an  inherent  error  detection  and  correction  capability  such  as 
analog-to-digital  converters,  direction  finding  arrays  and  radar  waveform  design.  To  use  the  RSNS,  it  is 
necessary  to  know  the  greatest  length  of  combined  sequences  without  ambiguities,  called  the  dynamic 
range  cM,  for  which  only  a  few  closed-form  expressions  currently  exist.  In  this  paper,  an  efficient 
algorithm  for  computing  cM  and  its  position  within  the  combined  set  of  sequences  is  presented  and  shown 
to  be  independent  of  the  size  of  the  moduli.  The  algorithm  is  used  to  generate  the  equations  for  several 
groups  of  additional  moduli  arrangements.  Closed  form  expressions  for  cM  are  conjectured  and  proved 
using  the  obtained  congruence  equations  that  define  the  ambiguity  locations. 

15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OF 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

li 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


This  article  has  been  accepted  for  publication  in  a  future  issue  of  this  journal,  but  has  not  been  fully  edited.  Content  may  change  prior  to  final  publication.  Citation  information:  DOI 

10.1 109/TIT. 2014. 2301820,  IEEE  Transactions  on  Information  Theory 

IEEE  TRANSACTIONS  ON  INFORMATION  THEORY,  VOL.  XX,  NO.  XX,  MONTH  YEAR  2 


expressions  from  the  congruence  equations. 

The  paper  is  structured  as  follows.  In  Section  II,  we  briefly 
review  the  RSNS  formulation.  In  Section  III,  RSNS  ambigui¬ 
ties  are  introduced.  In  Section  IV,  minimal  ambiguity  pairs  are 
defined  and  a  series  of  lemmas  and  theorems  are  presented  that 
provide  the  details  of  our  algorithm  and  the  efficiency  of  our 
solution.  An  TV  =  3  example  is  also  provided  to  demonstrate 
the  steps  of  the  algorithm.  In  Section  V,  the  algorithm’s 
run  time  complexity  is  compared  to  that  of  a  naive  search 
approach.  The  complexity  for  both  computations  is  derived 
demonstrating  the  algorithm  efficiency  and  its  independence 
from  the  size  of  the  TV  moduli.  In  Section  VI,  the  results  for 
several  groups  of  moduli  sets  are  used  to  generate  an  extended 
group  of  closed-form  expressions  for  M.  In  section  VII,  we 
provide  concluding  remarks  and  areas  of  further  research. 

Throughout  this  paper,  we  use  the  Vinogradov  symbols  2>, 
<  and  the  Landau  symbols  O,  o  with  their  usual  meanings. 
We  recall  that  /  <7,  g  3>  /  and  /  =  O(g)  are  all  equivalent 

and  mean  that  \f(x)\  <  c\g(x)\  holds  with  some  constant  c, 
for  x  sufficiently  large.  Also,  /  =  o(g)  if  linx^oo  =  0. 
For  a  positive  real  number  x  we  write  log  x  for  the  maximum 
between  2  and  the  natural  logarithm  of  x. 

II.  Robust  Symmetrical  Number  System 

The  RSNS  is  a  modular  based  number  system  consisting  of 
TV  >  2  integer  sequences  with  each  sequence  associated  with 
a  coprime  modulus,  rrii.  The  RSNS  is  based  on  the  following 
sequence: 

{zm}  =  [0,1,2,  1,  m,  2,1],  (1) 

To  form  the  TV— sequence  RSNS,  each  term  in  (1)  is  repeated 
TV  times  in  succession.  Therefore,  the  integers  within  one 
folding  period  of  a  sequence  are: 

{xm}  —  [0,  •  ■  ■ ,  0, 1, 1, . . . ,  1, . . . ,  to  -  1, . . . ,  m  -  1,  ^ 

TO,  .  .  .  ,  TO,  TO  —  1,  .  . .  ,  TO  —  1,  .  .  .  ,  1,  .  .  .  ,  1]. 

This  results  in  a  periodic  sequence  with  a  period  of  Pm  = 
2 Nm  [4],  [9].  Each  sequence  corresponding  to  nii  is  also 
shifted  left  (or  right)  by  Si  =  i  —  1  where  i  £  {1,2...,  TV} 
and  the  shift  values.  Si  =  {si,  S2, . . . ,  sjy},  form  a  complete 
residue  system  modulo  TV.  The  resulting  structure  of  the  TV 
sequences  ensures  that  two  successive  RSNS  vectors  (paired 
terms  from  all  N  sequences)  when  considered  together,  differ 
by  only  one  integer  resulting  in  an  acyclic  integer  Gray  code 
property  [4],  [15].  Although  the  RSNS  is  cyclic  and  has  integer 
Gray  code  properties,  it  differs  from  an  (?n,  TV)-Gray  code 
in  the  following  ways:  all  the  TV— tuples  in  an  RSNS  are 
not  distinct,  and  the  maximum  value  of  each  element  in  an 
TV— tuple  is  not  to,  but  rather  nii,  where  i  represents  the 
elements  of  the  TV— tuple  [0,  1,  ...,  TV  —  1]. 

Each  sequence  is  extended  periodically  with  period  2Nm 
as  xh+n2Nm  =  xh  where  n  £  {0,  ±1,  ±2, . . .}.  There¬ 
fore,  Xh  is  a  symmetrical  residue  of  ( h  +  n2Nm)  modulo 
2Nm.  An  integer,  h,  is  represented  by  a  vector,  X h  = 
[xi^h,  X2,h:  ■  ■  ■  >  XN,h]T,  °f  N  paired  integers  from  each  se¬ 
quence  at  h.  For  example,  a  left-shifted,  three-sequence  RSNS 
with  rrii  =  {3,4,5}  is  displayed  in  Table  I  and  Fig.  1.  The 


mi  =  3 


Fig.  1.  RSNS  structure  for  mi  =  {3,4,5}. 


integer,  h  =  5,  is  represented  by  the  vector,  Xh  =  [1,2,  2]T. 
Also  the  integer  Gray  code  property  is  evident. 

Since  the  integer  values  within  each  modulus  consists  of 
2Nm  integers,  the  symmetrical  residues  are  determined  by 
first  subtracting  an  integer  number  of  2Nm  integers  as 


rii  =  h  — 


h 

2  TV  rrii 


2  TV  rrii. 


(3) 


The  symmetrical  residue  Xh  is  then  calculated  as  [4],  [9] 

77-j 


Xh  = 


TV 

2Nrrii  +  TV  —  rii  +  Si  —  1 
TV 


Si  <m  <  Nim  +  Si  + 1 


,  Nrrii  +  Si  +  2  <  rii  <  2 TVtoj  + 


(4) 

The  TV-sequence  RSNS  is  periodic  with  a  fundamental  period 
of 

Pf  =  21 VAf,  (5) 


N 

where  M  =  {{  rat  is  the  dynamic  range  of  a  residue  number 

system  (RNS)~~ [3],  [4],  [  1 3] .  _ 

Closed-form  expressions  for  M  exists  for  only  few  specific 
cases.  In  [4],  a  closed-form  expression  for  a  TV  =  2  RSNS  is 
reported,  where 


M  = 


4mi  +  2to2  —  5, 

4toi  +  2rri2  —  2, 


when  to2  <  m\  +  2 
when  m2  >  mi  +  3 


(6) 


and  5  <  mi  <  m2 .  The  other  published  closed-form 
expression  for  M  is  when  TV  =  3  and  rrii  =  {m— 1,  m,  to+1} 
with  m  even  and  to  >  3  [4],  [13].  In  this  case, 

_ ^3  ig 

M  =  -to2  +  +  7.  (7) 


III.  RSNS  Ambiguities 

In  the  fundamental  period  of  an  TV-sequence  RSNS,  there 
are  three  ambiguity  types.  Type  0,  Type  1,  and  Type  2  that 
are  illustrated  in  Fig.  2.  Type  0  ambiguities  occur  periodically 
in  each  sequence.  Type  1  ambiguities  occur  across  the  folds 
of  each  waveform  in  each  sequence,  and  Type  2  ambiguities 
occur  due  to  the  TV  sequential  repeated  integers  within  each 
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TABLE  I 

N  =  3  RSNS  STRUCTURE  FOR  m,  =  {3,  4,  5}. 


h 
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16  .  .  . 
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Type  0 


Fig.  2.  Single  sequence  ambiguity  types  (N=3,  m=5). 


sequence.  Each  sequence  can  be  decimated  into  TV  subse¬ 
quences  where  each  subsequence  is  composed  of  values  from 
the  original  sequences  at  h  =  0  mod  TV,  h  =  1  mod  TV, . . . , 
and  h  =  (TV  —  1)  mod  TV.  By  examining  these  subsequences, 
the  Type  2  ambiguities  are  eliminated  leaving  only  Type  0 
and  Type  1  ambiguities  [4].  Subsequently,  the  complexity  of 
the  problem  of  determining  M  is  reduced.  Table  II  illustrates 
the  subsequence  stmcture  for  a  single  sequence  of  a  TV  =  3 
RSNS.  By  analyzing  the  parity  of  the  RSNS  in  Table  I,  it  can 
be  shown  that  the  parity  of  the  sequence  repeats  at  a  period 
of  2TV  [4],  [16]. 

The  ambiguity  locations  can  be  determined  by  solving  the 
congruence  equations  for  each  combination  of  Type  0  and 
Type  1  ambiguities  that  exist.  To  determine  M,  the  ambiguity 
locations  are  determined  by  solving  congruence  equations  for 
all  combinations  of  Type  1  and  Type  0  ambiguities,  and  the 
largest  span  of  unambiguous  values  is  identified,  which  is 
equal  to  M.  The  various  combinations  of  Type  1  and  Type  0 
ambiguities  are  referred  to  by  three  digit  Case  Numbers  where 
the  first  digit  refers  to  the  number  of  Type  1  ambiguities  that 
exist  and  ranges  from  zero  to  N.  The  second  digit  represents 
the  particular  assignment  of  Type  0  and  Type  1  ambiguities 
to  specific  sequences,  and  third  digit  of  the  case  number 
represents  the  subsequence  index  and  ranges  from  zero  to 
N—  1  [4].  For  example,  for  an  TV  =  3  RSNS,  in  Case  220,  the  2 
as  the  first  digit  specifies  that  there  are  two  Type  1  ambiguities 
(and  therefore  one  Type  0  ambiguity)  for  the  three  sequences. 
The  2  as  the  second  digit  signifies  that  the  particular  order  of 
the  ambiguities  is  the  second  largest  binary  value  (IOI2  =  5). 
The  0  in  the  third  digit  of  the  example  case  label  indicates 
that  the  ambiguities  are  computed  for  the  0th  subsequence 


only  (see  [4],  [16]  for  further  details).  Table  III  summarizes 
the  solution  to  the  congruence  equations  and  defines  the  CO  A 
for  all  possible  case  numbers  for  an  TV-sequence  RSNS  [4], 
[16]. 

In  Table  III,  the  value  of  hs  is  determined  by  solving  a  set  of 
congruence  equations  using  the  Chinese  Remainder  Theorem 
(CRT).  The  congruence  equations  are  generated  from  a  matrix 
that  is  formed  by  collecting  the  Type  1  symmetrical  residue 
numerators  into  a  matrix  form  with  each  column  index  corre¬ 
sponding  to  its  associated  modulus.  The  matrix  has  a  unique 
structure  where  the  first  column  is  [0,  — 1,  —2, . . . ,  —N  +  1]T 
and  the  subsequent  columns  are  generated  by  circular  shifting 
the  previous  column  up  and  incrementing  each  value  by  one 
[4].  The  shift  matrix  is  an  N  x  N  matrix  defined  as 


0 

0 

.  0 

-1 

-1 

-1  TV-1 

-2 

. 

to 

to 

-TV +  2 

-TV +  2 

2 

_— TV  +  1 

1 

1  1 

(8) 


The  congruence  equations  are  formed  based  on  the  sequences 
containing  the  Type  1  ambiguities  and  the  subsequences  that 
contains  the  ambiguities. 

IV.  Efficient  Algorithm  for  Computing  M 

We  develop  an  efficient  algorithm  to  efficiently  compute  M 
for  N  RSNS  integer  sequences  with  arbitrary  coprime  moduli, 
nii,  where  m,  >  2,  by  first  considering  all  the  minimal  pair 
ambiguity  locations  For  the  general  TV— sequence 

RSNS  case,  we  let  C  =  {(/ii,/i2)  |  0  <  Tii  <  T12  <  P/}, 
where  X ^  =  X^2.  A  pair  (fi.1,/12)  £  C  is  minimal  if  there 
does  not  exist  a  pair  (hi,  £  C  such  that  hi  <  h\  <  K2  < 

T12  and  if  the  shorter  sequence  length,  h-2  —  h-\  <  h?  — hi.  The 
largest  distance  between  consecutive  minimal  pairs  h2~hi  —  l 
is  the  M  and  hi  +  1  is  the  starting  ^position  of  M.  We  will 
also  demonstrate  in  Theorem  7  that  M  <  Pf.  It  then  follows 
that  M  <  M  =  Y[mi,  the  dynamic  range  of  the  RNS. 

In  [4],  it  was  demonstrated  that  the  distance  between 
ambiguous  vector  pairs  is  always  odd;  therefore,  we  define 
the  midpoint  between  the  ambiguous  vector  pairs  as  CO  A  = 
(/12  +  hi)  /2.  Given  two  minimal  pairs.  Pi  =  (/ii,Ti2)  £  C 
with  COAp1  and  P2  =  (h\ .  h'2)  £  C  with  COAp2  where 
COAp1  <  COAp2,  the  pairs  are  defined  as  consecutive  if  there 
does  not  exist  a  minimal  pair  P3  =  ( h" .  h2 )  £  C  with  COAp3 
such  that  COAp1  <  COAp3  <  COAp2.  Therefore,  if  (hi,  h2)  £ 
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TABLE  II 

Single  sequence  RSNS  structure  illustrating  three  subsequences  for  anjV  =  3  RSNS. 


rrii  =  3 

Xh 

0 

0 

0 

1 

1  1  2 

2 

2 

3 

3 

3 

2 

2 

2 

1 

1 

1 

0 

h  =  0  mod  3 

Xh 

0 

1 

2 

3 

2 

1 

0 

h  =  1  mod  3 

Xh 

0 

1 

2 

3 

2 

1 

h  =  2  mod  3 

Xh 

0 

1 

2 

3 

2 

1 

h  | 

0 

1 

2 

3 

4  5  6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

TABLE  III 

Summary  of  N-sequence  RSNS  ambiguity  expressions.  Ambiguities  exist  at  h  and  h  +  k.  Index  i  denotes  the  sequences  with  Type  1 

AMBIGUITIES  AND  THE  INDEX  j  DENOTES  THE  SEQUENCES  CONTAINING  TYPE  0  AMBIGUITIES. 


Case  Number 

h  is 

k  is  a  multiple  of 

COA 

010 

Anywhere  in  Pf 

N 

2N  n  rn* 

i=i 

None 

1X0 

k 

h  =  aNrrii  —  — 

2 

M 

2N - 

rrii 

aNrrii 

1XX 

k 

h  =  aNrrii  +  hs  —  — 

2 

7F? 

2N - 

rrii 

hs  +  aNrrii 

2.Y0  •  •  •  (JV  -  1)X0 

k 

h  =  aN  FT  rrii  —  — 

i  2 

2  N  n  m, 

3 

aN  rrii 

2XX  •  •  •  (JV  -  1)XX 

k 

h  —  aN  FT  rrii  +  hs  —  — 

i  2 

2JV  n  m, 

3 

hs  +  aN  rrii 

AGO 

k 

h  =  aN  n  mn  —  — 

71=1  2 

2N 

71=iV 

aN  n  mn 

71=1 

NIX 

n  =  IV  /c 

h  =  aN  mn  +  hs  —  — 

71=1  2 

2N 

n  —  N 

hs  +  aN  n  mn 

71=1 

C  and  (h\ ,  h'2)  £  C  are  consecutive  minimal  pairs,  then  the 
maximal  size,  M  =  (/12  —  1)  —  (hi  +  1)  +  1  =  h<i  —  hi  —  1, 
is  the  dynamic  range  of  the  RSNS.  Furthermore,  hi  +  1  is 
the  beginning  position  of  the  dynamic  range.  Since  M  is 
computed  using  consecutive  minimal  pairs  (ft.1,/12)  £  C  and 
(h\ ,  h'f)  £  C,  only  the  positions  of  the  minimal  pairs  that  can 
affect  the  length  of  M  are  required  to  be  computed  and  the 
rest  can  be  ignored. 

The  algorithm  for  computing  M  relies  on  a  number  of 
lemmas,  most  of  which  are  the  result  of  an  analysis  of  the 
locations  of  all  vector  ambiguities  provided  in  [4],  Table  III 
summarizes  the  TV— channel  RSNS  vector  ambiguity  locations. 
The  rows  in  Table  III  separate  the  locations  of  the  ambiguity 
pairs  into  seven  categories  based  on  the  type  of  ambiguity. 


A.  Theoretical  Basis  for  Algorithm 

The  basis  for  our  efficient  algorithm  for  determining  M 
is  presented  in  a  series  of  lemmas  and  theorems.  From  this 
theoretical  foundation,  the  steps  of  the  algorithm  are  developed 
and  presented  in  Section  IV-B. 

Lemma  1.  There  are  2 N  distinct  cases  of  repeated  ambiguity 
pairs,  each  with  a  different  ambiguity  length  and  COA  spac¬ 
ing.  All  but  one  of  the  2N  cases  have  N  subcases  that  have 
the  same  number  of  COAs  and  ambiguity  lengths  in  Pf,  but 
the  COA  for  each  of  the  subcases  is  shifted  by  a  particular 
value,  hSi. 

The  N  COA  shifts  (one  for  each  subcase)  are  computed  by 
solving  a  set  of  N  congruence  equations  using  the  CRT.  The 
subcase  where  hSQ  =  0  is  called  the  base  case  and  is  shown 
in  rows  2,  4,  and  6  in  Table  III. 


Proof.  See  the  ambiguity  analysis  discussion  in  Section  II 
of  [4],  □ 

Lemma  2.  Minimal  pairs  are  computed  using  the  first  multiple 
of  k  from  the  third  column  in  Table  III. 

Proof.  Any  vector  pair  computed  using  a  higher  multiple  of 
k  forms  a  vector  pair  that  encompasses  and  is  symmetric 
about  the  vector  pair  obtained  using  the  lower  multiple  of  k. 
Therefore,  any  vector  pair  computed  using  other  than  the  first 
multiple  of  k  is  not  minimal  [4],  □ 

Lemma  3.  For  every  ambiguity  pair  with  a  COA  at  h,  there 
is  an  ambiguity  pair  with  the  same  length  at  h  +  Pf  / 2. 

Proof.  Given  a  general  COA  for  any  case  at 


where  the  subscripts  i  are  the  indices  of  all  vector  elements 
with  Type  1  ambiguity,  there  is  also  a  COA  at  h+Pf  / 2  because 


(10) 
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where  a  and  b  are  any  integers,  j  corresponds  to  the  vector 
elements  with  Type  0  ambiguities,  and  i  corresponds  to  the 
vector  elements  with  Type  1  ambiguities.  The  result  is  that 
ambiguity  pairs  are  symmetric  about  Pf/2.  □ 

Lemma  4.  There  is  always  an  ambiguity  with  a  COA  at  h  =  0 
and  h  =  Pf  /2  with  a  length  of  2 N  +  1,  which  is  also  the 
ambiguity  with  the  smallest  length. 


For  the  case  N  =  3,  we  can  derive  an  exact  expression  for 
the  upper  bound  on  M,  which  we  do  in  the  next  lemma. 

Lemma  8.  In  the  case  of  3 -channel  RSNS  of  coprime  moduli 
TOi  <  m2  <  m3,  the  dynamic  range  is  upper  bounded  by 


— 


iV(mim,2  +  2mf)  —  1 
AT  (m3  +  2m\m2)  —  1 


if  mi  m2  >  m3 
if  m2m2  <  m3. 


03) 


Proof.  This  is  straightforward  via  inspection  of  row  6  in 
Table  III.  □ 

Lemma  5.  Using  Lemma  3  and  Lemma  4,  only  ambiguities 
from  —N  to  Pf  /2  +  N  need  to  be  considered  when  comput¬ 
ing  M. 

Proof.  Since  Lemma  3  showed  that  ambiguity  pairs  are 
symmetric  about  Pf/ 2,  and  M  is  computed  from  minimal 
ambiguity  pair  locations,  the  same  length  M  exists  from  h  =  0 
to  h  =  Pf/ 2  as  exists  from  h  =  Pf / 2  to  h  =  Pf.  □ 

Lemma  6.  [4]  The  dynamic  range  is  upper  bounded  by 


Proof.  We  need  to  minimize  the  expressions  (from  (11)): 

«i  =  mim.2  +  2m.3 
a2  =  mim3  +  2  m2 
c*3  =  m2m3  +  2m  1 

«4  =  mim.2m.3  +  2 

I 

0:5  =  1  +  2mi?7l2?77.3 

«6  =  m±  +  2m-2m.3 

«7  =  m2  +  2m-im.3 
a  8  =  m3  +  2m.im.2- 


M  <  B-rf  :=  N  min 

M  I  C{l....,n} 


ntj  +  2  mj 
iei  jet 


-1, 


(11) 


By  inspection,  it  is  easy  to  see  that  a4  <  a2  <  0.3  <  a 4  and 
a%  <  oit  <  otQ  <  0:5.  That  is,  we  need  to  only  compare  a\ 
and  as,  which  is  equivalent  to  comparing  m4m2  and  m3.  □ 


where  I  is  the  set  complement,  that  is,  I  =  {1, ...  ,n}\I. 

That  is,  each  of  the  rows  in  Table  III  produce  a  unique 
set  of  minimal  pairs  and  the  row  that  has  the  smallest  local 
M  (the  one  that  minimizes  (11))  provides  an  upper  bound 
on  M  for  the  RSNS.  Any  ambiguity  pair  that  has  a  length 
greater  than  ^  does  not  affect  M  and  can  be  ignored  (i.e., 
M  is  smaller  than  the  distance  between  the  minimal  pair  and 
therefore  cannot  be  the  vector  pair). 

As  an  example,  let  N  =  3,  m,:  =  {3, 4, 5}.  We  first  compute 
the  expressions  inside  the  minimum  of  (11),  for  each  of  the 
2"  subsets  /  C  {1,2,...,  n},  that  is, 

Pi  =  [1  +  2(3  -4-  5)]  =  121 
B2  =  [3  +  2(4  •  5)]  =  43 
B3  =  [ 4  +  2(3  •  5)]  =  34 
B4  =  [5  +  2(3  •  4)j  =  29 
P.5  =  [(3  •  4)  +  2(5)]  =  22 
P6  =  [(3  •  5)  +  2(4)]  =  23 
B7  =  [(4  •  5)  +  2(3)]  =  26 
B8  =  [(3-4-5)  +  2(l)]=65. 

By  using  Lemma  6, 

Bg  =  3  min  (Bt)  -  1  =  3  •  22  -  1  =  65. 

1V1  2 

Theorem  7.  Assuming  M  >  4,  we  can  take  as  upper  bound 
for  M 

B~  <  N\2V2M]  -  1.  (12) 

Moreover,  if  N  >  3,  the  dynamic  range  of  the  RSNS  is  always 
smaller  than  the  dynamic  range  of  the  RNS,  that  is, 

M  <  M. 

Proof.  See  Appendix.  □ 


For  instance,  if  N  =  3  and  m;  =  {3,4,  5},  since  mim2  > 
m3,  B^  =  3(3  •  4  +  2  •  5)  -  1  =  65. 

B.  Efficient  Algorithm  Steps 

Using  Lemmas  1  through  8,  the  efficient  algorithm  for 
computing  M  follows  the  steps  below: 

51.  Define  N  as  the  number  of  coprime  moduli  (m,)  i<,;</Y, 
M  =  Jl]=i  TOi’  al,3  fundamental  period  of  the  RSNS 

Pf  =  2  AIN. 

52.  Compute  the  upper  bound  Bff  for  the  dynamic  range, 
(12)  of  Theorem  7,  or  (13)  of  Lemma  8,  if  N  =  3. 

53.  Compute  the  number  of  ambiguity  cases  for  the  par¬ 
ticular  RSNS  using  Table  III  and  the  limits  imposed  by 
Lemma  1 .  Compute  the  minimal-pair  distance  for  all  am¬ 
biguity  pair  cases  using  multiplication  of  corresponding 
entries  in  the  matrix  of  size  N  x  2^,  which  is  the  matrix 
containing  as  columns  all  of  the  N  subcases  of  the  2N 
distinct  cases  of  repeating  ambiguity  pairs,  and  eliminate 
all  ambiguity  pair  cases  that  have  a  length  greater  than 
the  dynamic  range  upper  bound  (step  S2). 

54.  Compute  the  remaining  minimal  pair  ambiguity  loca¬ 
tions  (h\,h2)  using  Table  III  and  Lemmas  2  and  5. 

55.  Sort  the  matrix  of  minimal  pairs  (h\,  h2)  such  that  h2  is 
monotonically  increasing.  Vector  subtract  the  start  posi¬ 
tions  of  consecutive  minimal  pairs  (h\  ( p )  —  h4  (p  +  1)) 
and  remove  all  minimal  pairs  where  the  result  is  nega¬ 
tive.  The  remaining  minimal  pairs  are  the  only  consec¬ 
utive  minimal  pairs. 

56.  Compute  the  vector  of  distances  between  endpoints  of 
consecutive  minimal  pairs  (h2  (p  +  1)  —  h±  (p)  —  1).  The 
dynamic  range  M  is  the  largest  value  in  the  resulting 
vector. 
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C.  Example 

Consider  the  RSNS  from  Table  I  where  TV  =  3,  n%i  = 
{3,4,5}.  Table  III  provides  the  general  TV— sequence  ambi¬ 
guity  expressions  and  Table  I  of  reference  [4]  provides  the 
general  ambiguity  pair  equations  for  the  TV  =  3  case  (we  give 
a  particular  case  of  that  table  in  our  Table  IV,  for  the  moduli 
mi  =  {3,4,5}). 

51.  Define  TV  =  3,  M  =  fl^Li  m*  =  60,  and  Pf  =  360. 

52.  Compute  the  dynamic  range  upper  bound  by  using 
Theorem  7  (in  which  case,  [6\/120]  —  1  =  [65.7267]  — 
1  =  65),  or  (since  TV  =  3)  Lemma  8  (in  which  case, 
3(3  •  4  +  2  •  5)  -  1  =  65). 

53.  Table  IV  shows  all  possible  ambiguity  cases  (using 
Table  III,  or  Table  I  in  [4])  and  points  out  (above  the 
double  line)  the  rows  that  have  ambiguity  pairs  with  a 
length  (min  k)  greater  than  =  65,  which  can  be 
ignored  in  the  computation  of  M. 


TABLE  IV 

All  ambiguity  cases  (after  [4],  a  e  Z) 


Case 

Label 

Ambiguities  occur  at 
h  and  h  +  k,  where  h  is 

min  k 

Case  010 

Any  position 

360 

Case  110 

h-Case  110  =  a  ■  9  —  60 

120 

Case  111 

h-case  111  =  a  ■  9  —  59 

Case  112 

h-Case  112  =  a  ■  9  —  58 

Case  120 

h'Case  120  =  cl  •  12  —  45 

90 

Case  121 

h'Case  121  =  a  ■  12  —  44 

Case  122 

h'Case  122  —  CL  ■  12  —  46 

Case  130 

h'Case  130  =  a  ‘  15  —  36 

72 

Case  131 

he ase  131  =  CL  •  15  —  38 

Case  132 

he ase  132  =  CL  ■  15  —  37 

Case  210 

he  ase  210  =  cl  •  36  —  15 

30 

Case  211 

he  ase  211  =  a  •  36  +  hs  i  —  15 

Case  212 

he  ase  212  =  CL  •  36  +  hs  2  —  15 

Case  220 

he  ase  220  =  cl  ■  45—12 

24 

Case  221 

he  ase  221  =  a  •  45  +  hs  1  —  12 

Case  222 

he  ase  222  =  CL  •  45  +  hs  2  —  12 

Case  230 

h-Case  230  =  CL  ■  60  —  9 

18 

Case  231 

he  ase  231  =  CL  ■  60  +  hsi  —  9 

Case  232 

he  ase  232  =  CL  ■  60  +  hs  2  —  9 

Case  310 

he  ase  310  =  cl  •  180  —  3 

6 

Case  311 

hcase  311  =  CL  ■  180  hs  i  —  3 

Case  312 

he  ase  312  =  CL  •  180  +  hs  2  —  3 

The  base  cases  (all  cases  ending  in  zero  in  Table  IV) 
do  not  have  shifts  applied  to  the  CO  A  (i.e.,  hs  o  =  0). 
The  shifts  hs i  and  hs 2  in  Table  IV  are  computed, 
according  to  the  procedure  described  in  [4],  as  the 
least  positive  solutions  to  the  following  two  sets  of 
congruence  equations  (hs  1  =  73,  hs 2  =  119): 


ha  1  —  1 

=  0  mod  3 

fts  2  —  2 

=  0  mod  3 

3 

3 

hs  i  —  1 

=  0  mod  4 

fts  2  +  1 

=  0  mod  4 

3 

3 

hs  i  +  2 

=  0  mod  5 

hs  2  +  1 

=  0  mod  5. 

3 

3 

54.  Minimal  pair  ambiguity  locations  (hi,  /12)  are  computed 
using  Table  IV  for  h  =  —3  to  h  =  183.  All  minimal 
pairs  are  provided  in  Table  V. 

55.  The  consecutive  minimal  pairs  are  shown  in  Table  V. 

TABLE  V 

All  minimal  pairs,  and  all  consecutive  minimal  pairs  (shaded) 


hi 

COA 

h2 

hi 

COA 

h2 

-3 

0 

3 

68 

83 

98 

-14 

1 

16 

78 

90 

102 

-4 

11 

26 

93 

108 

123 

4 

13 

22 

94 

109 

124 

16 

28 

40 

106 

118 

130 

17 

29 

41 

116 

119 

122 

21 

36 

51 

111 

120 

129 

22 

37 

52 

124 

133 

142 

33 

45 

57 

123 

135 

147 

1  32 

47 

62  1 

129 

144 

159 

50 

59 

68 

130 

145 

160 

51 

60 

69 

140 

155 

170 

57 

72 

87 

151 

163 

175 

70 

73 

76 

152 

164 

176 

62 

74 

86 

170 

179 

188 

111 

180 

183 

S6.  The  consecutive  minimal  pairs  that  have  the  largest 
distance  between  them  are  displayed  in  bold  font  with  a 
shaded  background  in  Table  V.  The  result  is  an  M  =  43 
starting  at  h  =  79,  and  ending  at  ft,  =  121,  which  agrees 
with  the  results  in  [4],  [7], 

V.  Algorithm  Efficiency 

We  compute  the  complexity  of  the  computation  of  the  bound 
(11)  and  compare  the  efficiency  of  our  algorithm  with  the  naive 
search  algorithm  of  [3],  We  use  the  “prime”  big-oh  notation 
O'(-)  for  functions  in  both  M,  TV  (to  see  the  dependence  on 
TV),  and  the  big-oh  notation  O(-)  for  functions  in  M,  which  is 
the  relevant  parameter  (the  O-constant  will  be  dependent  on 
TV).  For  every  0  <  k  <  TV,  and  every  subset  I  of  cardinality  k , 
we  perform  k  —  1  multiplications  for  the  first  term  in  the  min¬ 
imum  computation  of  (11),  plus,  a  division  and  a  doubling  for 
the  second  term  (since  2Y[jeImj  =  2M/ J}igj reusing 
the  previous  computation).  Including  the  sorting  for  a  set  of 
cardinality  2N  with  complexity  0(N2n),  the  complexity  for 
the  bound  computation  of  Lemma  6  is 


=  O  (N2n  +  2N3n~1) 
=  O  (N3n)  , 


using  the  identity  J2k= oMfc)2^  =  nz(z  +  l)n_1.  Applying 
(12)  significantly  reduces  the  above  complexity  to  0(1),  as  a 
function  of  M  (at  the  expense  of  possibly  increasing  the  upper 
bound).  The  M  computation  process  described  in  Section  IV-B 
was  implemented  using  MATLAB.  A  search  of  current  re¬ 
search  did  not  reveal  any  existing  efficient  computational 
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Fig.  3.  Log-log  plot  of  M  vs.  run  time  using  the  improved  and  naive 
algorithms. 


algorithms  for  finding  and  comparing  all  pairs  of  N  x  1  vectors 
in  an  N  x  Pf  vector  space,  except  for  a  naive  search  algorithm 
used  in  [5].  The  naive  search  algorithm  starts  by  creating 
an  N  x  Pf  matrix  with  each  of  the  columns  consisting  of 
the  integer  values  within  each  RSNS  modulus  sequence,  as 
shown  in  Table  I.  A  double  nested  for-loop  then  determines 
the  beginning  position  of  each  ambiguity  hi,  which  are  then 
sorted.  A  second  double  nested  for-loop  is  then  used  to 
determine  the  end  position  /12  of  each  ambiguity  where  no 
other  ambiguities  are  enclosed.  The  maximum  length  is  then 
calculated  and  is  M  with  the  matrix  index  corresponding  to 
the  correct  beginning  and  end  positions. 

Now,  we  compare  the  time  complexities  (arithmetical  oper¬ 
ations,  and  comparisons)  of  both  algorithms  in  the  modulus  M 
(assuming  N  fixed).  The  naive  approach  uses  a  matrix  of  size 
N  x  (2 NM)  and  for  each  column,  it  checks  for  the  first  match 
(ambiguity)  in  the  remaining  columns  of  the  matrix,  so  it  uses 
N  comparisons  (for  each  components  of  every  vector)  plus 
an  addition  for  the  range  counter.  Therefore,  the  worst  case 
complexity  of  finding  ambiguities  and  the  distance  between 
them  is 

pf- 1 

<  ^2  N(Pf  -h)  =  NPj  -  NPf^  f  ~1’  (15) 

h= 0 

=  NPf(Pf  +  ^  =  M(2MN  +  1  )N2  =  0'{M2N ?$£) 

(since  the  RSNS  fundamental  period  is  Pf  =  2 NM).  We  then 
sort  the  obtained  list  of  size  <C  Pf  =  2 NM,  which  can  be 
done  in  O' (MN (log  M  +  log  N)),  resulting  in  a  total  time 
complexity  of 

Naive  Time  Complexity  =  0'(M2N3)  =  0(M2) 

for  the  dynamic  range  computation.  Now,  we  examine  the 
time  complexity  of  our  algorithm.  The  first  step  is  the  same 
for  both,  and  we  disregard  its  complexity  (as  it  is  quite  low. 


that  is,  O'(N),  compared  to  the  other  steps’  complexities). 
The  second  step  can  be  done  in  0(1)  using  our  Theorem  7, 
Lemma  8.  The  third  step  uses  a  sorted  bin  of  size  N  x2N  and 
performs  Or(N2  2N)  additions/multiplications  on  rows,  and 
0'(N 22N)  additions/multiplications  on  columns  (see  [4]  for 
further  details  on  this  step).  The  sorting  of  Step  3  is  done  in 
0'(N2n)  operations.  Step  4  as  well  as  Step  6  need  0'(2N) 
operations,  and  Step  5  needs  0'{N 2N),  for  a  total  worst  case 
time  complexity  of 

Improved  Time  Complexity  =  O' (N22N)  =  0(1). 

Remark  9.  The  main  advantage  of  our  algorithm  is  that  it 
removes  the  apparent  dependence  on  the  size  of  the  moduli 
in  the  number  of  operations  needed  to  compute  the  dynamic 
range. 

Fig.  3  shows  a  log-log  plot  of  M  versus  computation  time 
for  the  two  algorithms  using  hundreds  of  N— channel  moduli 
sets.  Each  “+”  represents  the  M  obtained  using  the  algorithm 
presented  in  this  paper,  and  has  a  corresponding  “o”  on  the 
same  horizontal  axis  (up  to  M  «  109),  which  is  the  M 
computed  using  the  naive  search  algorithm.  The  results  are 
displayed  where  the  two  computation  methods  produced  the 
same  results  for  M  (up  to  104s).  For  example,  with  N  =  4 
moduli  with  M  =  2  x  104,  the  naive  algorithm  takes  300  s 
to  produce  the  answer,  while  the  efficient  algorithm  described 
above  only  takes  0.02  s. 

VI.  Further  Closed-Form  Results  for  M 

In  this  section,  closed-form  expressions  are  developed  for 
M  for  several  groups  of  moduli  sets  by  curve  fitting  data 
generated  using  the  efficient  algorithm  described  in  Section  IV 
and  then  verifying  that  the  closed-form  expressions  satisfy 
the  ambiguity  equations  of  Table  III.  Curve  fitting  the  data 
obtained  from  the  algorithm  was  chosen  as  a  method  of 
generating  the  closed-form  expressions  for  M  because  it 
leveraged  the  efficiency  of  our  algorithm  to  generate  data 
for  a  large  groups  of  moduli  sets  and  curve  fitting  enabled 
determining  patterns  of  repeated  Start  and  Stop  Case  number 
combinations  and  visually  determining  the  periodic  nature  of 
discontinuities  in  the  plotted  data  and  resolving  them.  This 
method  also  proved  to  be  efficient  compared  to  attempting  to 
derive  closed-form  expressions  for  M  by  iteratively  solving 
the  equations  of  Table  III  for  different  moduli  sets. 

A  sample  of  the  data  generated  using  the  algorithm  is 
presented  in  Table  VI.  The  data  was  analyzed  to  determine 
which  moduli  sets  have  the  same  Start  and  Stop  Case  numbers, 
and  the  value  of  M  was  plotted  against  a  variable  m  that 
is  linearly  related  to  the  moduli.  Curve-fitting  of  the  plotted 
data  using  MATLAB’s  curve  fitting  tool  box  was  conducted 
to  generate  exact  closed-form  expressions  for  M.  The  closed- 
form  expressions  were  then  verified  to  be  accurate  by  deriving 
them  from  the  applicable  equations  listed  in  Table  III.  Several 
groupings  of  moduli  for  N  =  3  and  N  =  4  were  examined 
increasing  the  number  of  closed-form  expressions  for  M. 

For  the  N  =  3  case,  sequential  coprime  odd  moduli  of  the 
form  rrii  =  {m  —  1,  m  +  1,  to  +  3}  with  m  even  and  m  >  3 
were  examined.  Also,  the  case  of  two  odd  moduli  and  one  even 
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modulus  were  examined  where  m*  =  {to,  to  +  1,  m  +  3}  and 
{to  —  3,  m  —  1,  to}  with  m  even.  The  final  case  examined  was 
where  the  moduli  consisted  of  every  other  odd  number,  that 
is  TOj  =  {to,  to  +  4,  to  +  8}  where  to  is  odd  and  to  >  3. 
The  new  closed-form  expressions  derived  from  curve  fitting 
the  data  are  presented  in  Table  VIII. 

To  demonstrate  the  method  used  to  verify  the  closed-form 
expressions  generated  from  curve  fitting,  we  examine  the  case 
where  N  =  3  and  rrii  =  {to  —  1,  to  +  1,  to  +  3}  with  to  even 
by  generating  the  closed-form  expressions  from  the  equations 
of  Table  III.  The  data  was  examined,  and  it  was  determined 
that  two  distinct  sets  of  case  numbers  are  associated  with  the 
beginning  and  ending  positions  of  M.  When  m  =  0  mod  4, 
the  case  number  associated  with  the  beginning  position.  Start 
Case  of  M  is  Case  211,  and  the  case  number  associated  with 
the  ending  position,  Stop  Case  of  M  is  Case  220.  The  value 
of  to  was  plotted  against  the  value  of  M,  and  the  data  was 
curve  fitted  to  a  quadratic  polynomial  using  MATLAB’s  curve 
fitting  toolbox.  From  the  curve  fit  data, 

M  =  jj to 2  +  -^to  +  7  (17) 

where  to  =  0  mod  4.  A  sample  of  the  data  used  to  derive 
(17)  is  presented  in  Table  VI.  The  data  examined  and  its 
corresponding  curve  fit  is  illustrated  in  Fig.  4a.  The  values  of 
M  derived  from  the  algorithm  are  equal  to  the  values  resulting 
from  (17). 

Now,  we  will  verify  that  the  curve  fit  solution  satisfies  the 
ambiguity  equations  of  Table  III.  For  the  Start  Case, 

h2n  =  a  (3?tiito2)  +  hS211  —  3to3.  (18) 


The  congmence  equations  that  are  generated  from  the  shift 
matrix,  (8), 

^ s 211 — -  =  0  mod  TOi 

h  -i  <19) 

S21\ - =  0  mod  TO2, 


are  solved  resulting  in  hS211  =  1.  The  value  of  a  in  (18)  was 
derived  by  solving  (18)  for  a  using  the  values  of  the  beginning 
position  (hi  +  1)  of  M  and  to.  The  results  were  curve  fitted, 
resulting  in  a  =  0.75m +  2.  By  substituting  the  expression  for 
a  into  (17),  we  obtain 


9  21 

/1211  =  -to3  +  6  to2 - to  —  14. 

4  4 


For  the  Stop  Case, 


ft-220  =  a(3TOi?n3)  +  3  to2. 


(20) 


(21) 


It  was  determined  that  a  =  0.75 to  +  1  and 

h220  =  +  y to 2  +  ^to  -  6  (22) 

By  solving  for  M  =  /1220  —  /1211  —  1,  we  obtain  (17). 

When  to  =  2  mod  4,  the  Start  Case  is  Case  211,  and  the 
Stop  Case  is  Case  231  as  shown  in  Table  VII.  The  data  was 
curve  fitted  resulting  in 

—  3  n  15 

M  =  -TO2  +  yTO  +  5.  (23) 


(a) 


(b) 


Fig.  4.  Curve  fitting  results  for  M  when  (a)  m,  =  {m  —  l,m+l,Tn  +  3} 
where  m  =  0  mod  4,  (b)  m;  =  {m  —  1,  m  +  1,  m  +  3}  and  m  =  2  mod 

4. 


Fig.  4b  displays  the  data  and  closed-form  solution  generated 
from  curve  fitting.  Using  the  same  approach,  (23)  is  verified 
by  deriving  it  from  the  equations  in  Table  III.  For  the  Stop 
Case, 

ft-231  =  a  (3 to2to3)  +  h.S231  +  3m  1 .  (24) 


The  congruence  equations  generated  from  (8) 

^ s 231 -  =  0  mod  to. 2 

o 

hs231  +  2  =  0  mod  to3  . 

3 

are  solved  using  the  CRT,  resulting  in 

,  3  2  15 

^231  =  +  ym  +  7. 


(25) 


(26) 


The  value  of  a  in  (24)  was  determined  to  be  a  =  0.5 to  —  1. 
These  expressions  were  then  substituted  into  (24)  to  determine 
that 

3  9 

ft23i  =  7; rn 3  +  -to2  +  3m  -  5.  (27) 

After  solving  for  M  =  ft.23i  —  /1211  —  1,  we  obtain  (23)  which 
verifies  the  result.  The  closed-form  expressions  for  the  other 
moduli  sets  examined  were  verified  in  a  similar  manner. 

Several  groups  of  N  =  4  RSNS  moduli  sets  were  also 
examined,  and  closed-form  expressions  to  M  were  produced 
by  curve  fitting  the  data  generated  by  the  efficient  algorithm. 
The  moduli  sets  examined  were 


TOi  =  {TO  —  1,  TO,  TO  +  2,  TO  +  4}  , 
TOi  =  {m,  TO  +  1,  TO  +  2,  TO  +  4}  , 
TOi  =  {TO,  TO  +  2,  TO  +  3,  TO  +  4}  , 
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TABLE  VI 

Sample  of  data  used  in  curve  fitting  for  sequential  odd  coprime  moduli,  m  =  0  mod  4 


m  —  1 

m  +  1 

m  +  3 

h±  +  1 

Start  Case 

h.2  —  1 

Stop  Case 

M 

7 

9 

11 

1481 

211 

1643 

220 

163 

11 

13 

15 

4676 

211 

4988 

220 

313 

15 

17 

19 

10655 

211 

11165 

220 

511 

19 

21 

23 

20282 

211 

21038 

220 

757 

23 

25 

27 

34421 

211 

35471 

220 

1051 

27 

29 

31 

53936 

211 

55328 

220 

1393 

31 

33 

35 

79691 

211 

81473 

220 

1783 

35 

37 

39 

112550 

211 

114770 

220 

2221 

39 

41 

43 

153377 

211 

156083 

220 

2707 

43 

45 

47 

203036 

211 

206276 

220 

3241 

47 

49 

51 

262391 

211 

266213 

220 

3823 

TABLE  VII 

Sample  of  data  used  in  curve  fitting  for  sequential  odd  coprime  moduli,  with  m  =  2  mod  4. 


m  —  1 

m  +  1 

m  +  3 

h±  +  1 

Start  Case 

h2  —  1 

Stop  Case 

M 

5 

7 

9 

395 

211 

498 

231 

104 

9 

11 

13 

1745 

211 

1974 

231 

230 

13 

15 

17 

4631 

211 

5034 

231 

404 

17 

19 

21 

9629 

211 

10254 

231 

626 

21 

23 

25 

17315 

211 

18210 

231 

896 

25 

27 

29 

28265 

211 

29478 

231 

1214 

29 

31 

33 

43055 

211 

44634 

231 

1580 

33 

35 

37 

62261 

211 

64254 

231 

1994 

37 

39 

41 

86459 

211 

88914 

231 

2456 

41 

43 

45 

116225 

211 

119190 

231 

2966 

45 

47 

49 

152135 

211 

155658 

231 

3524 

TABLE  VIII  _ 

New  closed-form  expressions  for  M  for  N  =  3  RSNS. 


mi 

M 

m 

{m  —  1,  m  +  1,  ra  +  3} 

im2  +  yin  +  7 

m  =  0  mod  4 

|  m2  +  ^-m  +  5 

m  =  2  mod  4 

{m,  m  +  1,  m  +  3} 

§  m2  +  +  6 

m  =  2  mod  4  and  m  >  14 

{m  —  3,  m  +  1,  m} 

|  m2  +  |  m 

m  is  even  and  m  ^  6k  where  k  —  1,2,... 

f  m2  +  Mm  +  48 

m  =  1  mod  8 

{ra,  m  +  4,  ra  +  8} 

§  m2  +  +  35 

m  =  3  mod  8 

§m2  +  ^  +  34 

m  =  5  mod  8 

f  m2  +  +  45 

m  =  7  mod  8 

Tfii  =  {to,  to  +  2,  m  +  4,  m  +  5}  , 

and 

TOj  =  {?7l,  TO  +  2,  TO  +  4,  TO  +  6}  , 

where  to  is  odd.  The  closed-form  expressions  for  M  are 
presented  in  Table  IX  and  were  verified  in  the  same  manner 
as  the  N  =  3  cases. 

VII.  Concluding  Remarks 

This  paper  presents  an  algorithm  to  compute  M  and  de¬ 
termine  its  location  in  the  RSNS  sequence.  The  algorithm 
reduces  the  computation  time  by  several  orders  of  magnitude 
compared  to  a  previously  reported  naive  search  algorithm.  In 
addition,  we  demonstrate  that  our  efficient  algorithm  removes 
the  apparent  dependence  on  the  size  of  the  moduli  in  the 


number  of  operations  needed  to  compute  M.  Linear  inter¬ 
polating  the  data  in  Fig.  3,  it  would  take  the  naive  search 
algorithm  more  than  32  years  to  find  M  for  the  same  N  =  8 
sequence  (of  approximately  28  bit  moduli)  that  the  efficient 
algorithm  computed  in  approximately  30  seconds.  Moreover, 
the  efficient  algorithm  uses  far  less  memory  than  the  naive 
search  algorithm;  therefore,  it  can  be  used  to  determine  the 
N  sequence  M  and  position  for  moduli  sets  with  much  larger 
fundamental  periods. 

The  algorithm  was  also  used  to  generate  data  sets  from 
which  curve  fitting  produced  additional  closed-form  expres¬ 
sions  for  M  increasing  the  groups  of  moduli  sets  for  which 
analytical  expressions  for  M  exist.  The  new  closed-form  ex¬ 
pressions  for  M  significantly  increase  the  number  of  analytical 
expressions  that  are  known  for  M  greatly  increasing  the  utility 
of  the  RSNS  to  solving  many  types  of  engineering  problems. 
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TABLE  IX 

\I\Y  CLOSED-FORM  EXPRESSIONS  FOR  M  FOR  N  =  4  RSNS. 


mi 

M 

m 

{m  —  1,  m,  m  +  2,  m  +  4} 

10m2  +  6m  +  20 

m  =  3  mod  6  or  m  =  5  mod  6  with  m  >  15  and  m  7^  {29,  33} 

{m,  m  +  1,  m  +  2,  m  +  4} 

10m2  +  22m  +  20 

m  =  1  mod  6  or  m  =  3  mod  6  and  m  >  15 

{m,  m  +  2,  m  +  3,  m  +  4} 

10m2  +  30m  —  6 
10m2  +  30m  —  8 

m  =  {25  +  12fc,  29  +  12fe,  31}  where  k  =  0,  1,  2,  •  •  •  , 
m  =  {43  +  12fc,  47  +  12fc}  and  k  =  0, 1,  2,  •  •  • 

{m,  m  +  2,  m  +  4,  m  +  5} 

lOm-2  +  54m  +  20 

m  >  39,  gcd(m,  5)  =  1,  and  m  is  odd 

{m,  m  +  2,  m  +  4,  m  +  6} 

10m2  +  38m  +  56 

m  >  13  and  gcd(m,  3)  =  1 

Further  research  opportunities  are  investigating  additional 
moduli  sets  where  closed-form  expressions  exist,  developing 
a  general  closed-form  expression  for  M,  and  formulating 
an  overarching  theoretical  framework  that  relates  the  various 
symmetrical  number  systems  to  each  other. 

Appendix 

Proof  of  Theorem  7 

Proof.  Setting  x  =  mi<  we  see  that  the  expression  that 

must  be  minimized  in  the  right-hand  side  of  (11)  is  in  fact  the 
function 

/( x)  =  x  H - ,  x  >  1. 

x 

By  examining  the  derivative  of  f(x),  a  simple  calculus 
analysis  reveals  that  the  function  has  a  global  minimum  at 
x  =  \J2M,  namely  2\j2M,  and  the  first  inequality  is  shown. 

The  work  in  [4]  contains  the  proof  of  M  <  M  for  N  =  3. 
As  examples,  if  N  =  3,  the  smallest  size  coprime  moduli  sets 
are  listed  in  Table  X  in  lexicographical  order. 

_TABLE  X 

Examples  of  MRsns  and  M  for  N  =  3 


mi 

MrsnS 

M 

{2,3,5} 

28 

30 

{2,3,7} 

35 

42 

{2,3,11} 

46 

66 

{3,4,5} 

43 

60 

Now,  assume  that  N  >  4.  We  need  to  prove  that 
N\2y/2M~\  <  M,  for  N  >  4.  Starting  with  the  simple 
inequality 

N\2V2M 1  <  N(2y/2M  +  1), 

it  is  sufficient  to  show  that  N(2\/2M  +  1)  <  M,  which  is 
equivalent  to  N\/SM  <  M  —  N,  and  8N2M  <  M 2  +  N2  — 
2MN,  that  is,  M2  -  2N{4N  +  1  )M  +  N2  >  0.  Looking  at 
the  previous  inequality  as  the  sign  of  a  concave-up  parabola 
in  M,  we  see  that  the  inequality  is  true,  as  long  as 

M  >  4 N2  +  N  +  2N\/ 4N2  +  2 N.  (28) 

If  N  =  4,  then  M  >  210,  and  the  right-hand  side  of  (28) 
is  ~  209.881;  if  N  =  5,  then  M  >  2310,  and  the  right-hand 
side  of  (28)  is  ~  299.88.  Thus,  the  inequality  (28)  is  true  for 

4  <  N  <  5. 


Next,  it  is  observed  that  for  any  N  >  3  (use  the  fact  that 
the  moduli  are  coprime):  if  N  =  3,  then 

r  2  •  3  •  5  =  30,  N  =  3 
M  >  <  2  •  3  •  5  •  7  =  210,  N  =  4 
[efc. 

M  >  2  ■  3  •  5  =  30;  if  TV  =  4,  then  M  >  2  •  3  •  5  •  7  =  210, 
etc.  For  arbitrary  N,  an  easy  inductive  procedure  reveals  that 

M  >  PN #, 

where  P^v#  =  TTfeLiPfe  the  primorial  function,  and  pk 
is  the  fcth  prime.  It  is  well-known  (and  easily  derivable  by 
using  the  prime  number  theorem  [17])  that  P^v#  =  exp[(l  + 
o(l))  nlog n]. 

Assume  N  >  6.  It  is  immediate  that 

Pjv#  >  {N  +  1)!. 

(For  the  interested  reader,  a  better  asymptotic  estimate  Pn#  = 
exp[(l  +  o(l))  AT  log  IV]  is  well-known  [17].)  We  now  show 
that  for  N  >  6,  the  right  hand  side  of  the  above  inequality 
satisfies 

(N  +  1)!  >  4 N2  +  N  +  2N \/ 4N2  +  2 N,  (29) 

which  will  imply  our  claim. 

We  will  prove  (29)  by  induction  on  N.  If  N  =  6,  then 

(6  +  1)!  -  (4  •  62  +  6  +  2  •  6\/ 4  •  62  +  2  •  e)  >  4740, 

and  so,  the  inequality  (29)  is  true  in  this  case.  Assume  that 
the  inequality  is  true  for  N  and  we  show  it  for  N  +  1,  that 
is,  we  start  from  (29)  and  multiply  by  N  +  2  on  both  sides, 
to  get 

(N  +  2)!  >  {N  +  2)  •  [4N2  +N  +  2N\/4N2  +  2iv)  . 

It  will  be  sufficient  to  show  that 

(N  +  2)  [AN2  +  N  +  2N \J AN2  +  2N^  >  A(N  +  l)2 
+{N  +  1)  +  2  (N  +  l)\f  A(N  +  l)2  +  2(AT  +  1). 

Since 

{N+2)(4N2+N)-4(N+1)2-{N+1)  =  4N3+5N2-7N-5 

is  increasing  and  greater  than  0,  for  N  >  6,  the  previous 
inequality  will  follow  if 

{N +2)2N \/ AN2  +  2N  >  2(iV+l) y/4 (N  +  l)2  +  2 (N  +  1), 
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which  by  squaring  both  sides  transforms  into 

-24  -  88 TV  -  1207V2  -  407V3  +  807V4  +  727V5  +  167V6  >  0, 

which  is  certainly  true  for  TV  >  6.  This  concludes  the  proof 
of  Theorem  7.  □ 
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